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In the study of trapped two-component Bose gases, a widely used dynamical protocol is to start from the 
ground state of a one-component condensate and then switch half the atoms into another hyperfine state. The 
slightly different intra-component and inter-component interactions can then lead to highly nontrivial dynamics. 
We study and classify the possible subsequent dynamics, over a wide variety of parameters spanned by the trap 
strength and by the inter- to intra-component interaction ratio. A stability analysis suited to the trapped situation 
provides us with a framework to explain the various types of dynamics in different regimes. 
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I. INTRODUCTION 

Two-component Bose-Einstein condensates (BECs) are in- 
creasingly appreciated as a rich and versatile source of intri- 
cate non-equilibrium pattern dynamics phenomena. In addi- 
tion to experimental observations (il-O, pattern dynamics in 
two-component BECs has attracted significant theoretical in- 
terest (see, e.g., |[T4l - [28ll and citations in (l^ ). 

In a number of two-component BEC experiments reported 
over more than a decade, a standard technique has been to start 
from the equilibrium state of a single-component BEC, e.g., 
populating a single hyperfine state of ^^Rb, and then using 
a tt/ 2_pulse to switch half the atoms to a different hyperfine 
state • This results in a binary condensate where the two 
intra- species interactions (gu and ^22) and one inter- species 
interaction (gi2) are all slightly different from each other, but 
the starting state is the ground state determined by gu alone. 
Since it has been realized several times in several different 
laboratory setups, this is a paradigm non-equilibrium initial 
state for binary condensate dynamics. A thorough and gen- 
eral analysis of the dynamics subsequent to such a 7r/2 pulse 
is thus clearly important. In this article we present such an 
analysis, clarifying the combined role of the inter-species in- 
teraction (^12) and the strength A of the trapping potential. We 
provide a stability analysis mapping out regions of the X-gi2 
parameter space hosting different types of dynamics. Since it 
is now routine to monitor real-time dynamics in such experi- 
ments (e.g. fS]), we also directly analyze the real-time evolu- 
tion after a 7r/2 pulse. 

It is widely known that the ground state of a uniform two- 
species BEC is phase separated or miscible depending on 
whether or not the inter- species repulsion dominates over the 
self-repulsions of the two species, i.e., if 



^11^22 < (^12)' 



(1) 



then the ground state is phase separated |15]. This criterion 
is also a key ingredient in understanding dynamical features 
such as pattern dynamics in the density difference between the 



two species — such "spin patterns" emerge when the phase 
separation condition is satisfied. This can be understood as 
the onset of a modulation instability |[T6l - [T8l] . identified by the 
appearance of an unstable mode in the excitation spectrum 
around a reference stationary state. For a homogeneous situa- 
tion, linear stability analysis shows that modulation instability 
sets in when the condition of Eq. ([T]) is satisfied |[T6l - [T8l] . 

The situation is different in the presence of a trapping po- 
tential. Phase separation in the ground state, as well as the ap- 
pearance of modulation instability when starting from a mixed 
state, now requires larger inter- species repulsion fT9ll . 
This suggests that the region of parameter space where pat- 
tern dynamics occurs also depends on the trap. A trap is al- 
most always present in cold-atom experiments, and it is easy 
to imagine experiments where the trapping potential is not ex- 
tremely shallow but varies between tight and shallow limits. 
It is thus necessary to examine the relevance of Eq. ([T]) for 
trapped binary BECs. To this end, we explore different trap 
strengths spanning several orders of magnitude, and identify 
the appropriate extensions of Eq. ([T]) for the type of spin dy- 
namics resulting from the 7r/2 protocol described above. 

We focus on the effects of two parameters. First, we study 
effects of changing cross-species interaction ^12, thus gener- 
alizing Eq. ([T]) for trapped situations. Second, we explore the 
role of the relative strength of the trap with respect to the inter- 
actions. Our analysis, performed for a one-dimensional (ID) 
geometry, sheds light on the situation where gu and ^22 are 
close but unequal: (a) the stability analysis is performed for 
911 = 922 and their difference serves only to select appropri- 
ate instability modes; (b) the simulations are performed with 

922/911 = l-Ol- 

In Section ini we introduce the formalism and geometry. In 
Section [Till we show results from a linear stability analysis 
for a sequence of trap strengths, and identify and analyze rel- 
evant modulation instabilities. Through an analysis of unsta- 
ble modes, we present a classification of the parameter space 
into dynamically distinct regions, in relation to the prototyp- 
ical initial state explained above. This may be regarded as a 
dynamical "phase diagram". A remarkable aspect is that the 
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"phase transition" line most relevant to spin pattern dynam- 
ics does not arise from the first modulation instability (studied 
in Ref. |[T4ll ). This first instability mode is antisymmetric in 
space, and as a result is not naturally excited in a symmet- 
ric trap with symmetric initial conditions. Complex dynamics 
(not due to collective modes but rather due to modulation in- 
stability) is generated only when the first spatially symmetric 
mode becomes unstable, which occurs at a higher value of gi2. 

In Section [Ivl we provide a relatively detailed account of 
the time evolution. For each trap strength A, for values of 
gi2 not much larger than gu, we observe simple collective 
modes. Above a threshold value of ^12, the oscillation ampli- 
tude becomes sharply stronger, and at the same time the mo- 
tion becomes notably aperiodic, signaling that the dynamics 
is more complex than a combination of a few modes. Dy- 
namical spin patterns start appearing at this stage and become 
more pronounced as gi2 is increased further. The threshold 
value at which the dynamics changes sharply corresponds to 
the second modulation instability line rather than the first, as 
we demonstrate through careful choice of parameters in each 
region of the phase diagram derived from stability analysis. 

Some further connections between the stability analysis and 
dynamical features, relating to the length scale of generated 
patterns, appear in Section jV] In the concluding Section [Vll 
we place our results in context and point out open questions. 
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FIG. 1. (Color online.) Results from stability analysis. Squared 
eigenvalues cj^ of the stability matrix M are plotted against gi2, for a 
tight trap (top) and for a shallow trap (bottom left). The arrows show 
the values of gi2 for onset of the two instabilities, namely gi2 (onset 
of spatially antisymmetric modulation instability) and ^12 (onset of 
spatially symmetric instability). Typical eigenvectors corresponding 
to these two modes are shown in the panels on lower right. 



II. GEOMETRY AND FORMALISM 

The relevant time-resolved experiments have been per- 
formed in both quasi- ID geometries (highly elongated traps 
with strong radial trapping) and in a 3D BEC of cylindri- 
cal symmetry with the radial variable playing analogous role 
as the ID coordinate 10, [H]. Since the basic phenomena are 
very similar, we expect the same theoretical framework to de- 
scribe the essential features of each case. For definiteness, in 
this work we show results for ID geometry. We expect the 
general picture emerging from this work to be qualitatively 
true also for other geometries exhibiting the same type of spin 
dynamics. 

We describe the dynamics in the mean field framework at 
zero t emp erature, i.e., by two coupled Gross-Pitaevskii equa- 
tions 1I3QM32II: 

idt^i = - + |AV +^n|V^ip +^i2|V^2|')v^i ,(2) 

idt^2 =(^-^dl^ ^AV +^12|V^l|' +^22|V^2|')V^2 .(3) 

Condensate wave functions i/ji (x, t) and ^^2(2^, t) are normal- 
ized to unity, and A is the strength of the harmonic trap. Fac- 
tors of particle number and radial trapping frequency are ab- 
sorbed as appropriate into the effective ID interaction param- 
eters gij 1 6, 32, 33]. We consider purely non-dissipative dy- 
namics, i.e., we do not attempt to model experimental loss 
rates with a phenomenological dissipative term as done in, 
e.g.,Refs. 

The equations above are in dimensionless form because we 
measure lengths in units of trap oscillator length and time in 



units of inverse trapping frequency, for a hypothetical trap of 
unit strength (A = 1). The scale for trap strengths is itself 
fixed by imposing gn = 1. With this convention, small val- 
ues of A correspond to a BEC in the Thomas-Fermi limit. For 
comparison, we note that the parameters of the experiment of 
Ref. [6] corresponds to A of order 10~^ in these units. Of 
course one can switch between different units via the transfor- 
mation: X x/l, t t/P, A XP, g gU and?/^ ipVl, 
where / is an appropriately chosen scale. 

The initial state after a 7r/2 pulse involves both components 
occupying the ground state of a single-component system of 
interaction 2^ii, because the atoms were all in the first hyper- 
fine state before the pulse. We model this initial situation as a 
two-component BEC with gu = ^22 = 9i2- The 7r/2 pulse 
may then be regarded as a sudden change (a quantum quench 
1I29I] ) of the interaction parameters ^22 and gi2. 

We use gii = ^22 for the stability analysis of Section 
imi For the explicit time evolution reported in Section [IVl 
we use gii and ^22 values close but unequal: gu = 1, 
^22 = 1-01. This choice of close values is convenient for illus- 
trating the structure of the phase diagram, especially for shal- 
low traps. In rubidium experiments the difference between 
gu and ^22 is somewhat larger (in the common case using 
^'^Rb hyperfine states |1) = |F = l,mF = -1) and |2) = 
|F = 2,mi? = l)); however our insights should be relevant 
to a broad regime of possible experiments. A full exploration 
of the regime of arbitrary differences (gu — ^22) remains an 
open task, beyond the scope of the present manuscript. 

Numerical simulations presented in Section |IV] were per- 
formed using a semi-implicit Crank-Nicolson method 11341 

m. 
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III. STABILITY ANALYSIS AND DYNAMICAL "PHASE 
DIAGRAM" 

We provide in this section a stability analysis for gn = 
g22 that maps out the regions of X-gi2 parameter space which 
support pattern formation instabilities. 

Ideally, one might like to perform a stability analysis 
around the initial state. However, in contrast to the homo- 
geneous case 1 16], we are faced with the situation that the 
initial state is not a stationary state of the final Hamiltonian. 
The choice of reference state is therefore a somewhat subtle 
aspect of the present analysis. 

We use as reference state V^o(^) the lowest-energy spa- 
tially symmetric stationary state of the case gu = ^22, with 
parameter gi2 set to its final value. (For large ^12, this is 
not the ground state for these parameters, which is phase- 
separated.) This reference state has the advantage of looking 
relatively similar to our actual initial state (two components 
totally overlapping in space), and of being a stationary state 
of the Hamiltonian for which we analyze linear stability. Our 
reference state can be regarded as placing both components in 
the single-component ground state for interaction gu -\- gi2- 
We are not aware of a suitable stationary state even more sim- 
ilar to the actual initial state. We will see that our stability 
analysis around this reference state will predict remarkably 
well the main observed time-evolution features described in 
Section HVl 

Note that it is not natural to use ^117^ ^22, because station- 
ary states for such a case typically do not overlap completely 
in space. Instead, in our approach the difference between gu 
and ^22 will play the important role of selecting certain in- 
stability modes over others. For this reason, inferences from 
the present analysis apply only to small relative differences 
between gu and ^22- 

We linearize Eqs. © and © around the reference station- 
ary state V^o(^)- 

^l)2{x,t) = [V^o(^) + (5V^2(^,^)] exp(-z/it), (4) 

where ji is the chemical potential corresponding to the ref- 
erence state. By keeping only terms of the first order in 
5il)i{x^t) and 6ilj2{x,t), we obtain a system of linear equa- 
tions which can be cast in the form: 

Here is a matrix differential operator which, upon dis- 
cretization or upon expansion in a set of orthogonal functions, 
becomes the so-called stability matrix (e.g., 122, 24]). We an- 
alyze below the eigenmodes of the stability matrix, which we 
have obtained by numerically calculating the reference sta- 
tionary state V^o(^) and expanding in the basis of harmonic 
trap (non-interacting) eigenstates. 

Since we use gu = ^22 for the stability analysis, eigen- 
modes will have well-defined "species parity", i.e. will all 
be either even [Silji{x^t) = Silj2{x^t)] or odd [6ilji{x,t) = 
—Sijj2{xj t)] with respect to the interchange of species. Even 
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FIG. 2. (Color online.) The dynamical "phase diagram" showing 
the critical values of gi2 for the onset of the two types of modulation 
instability versus the trap strength A. The instability lines are shown 
as straight lines joining numerically determined values. The oscilla- 
tion schematics on the right (and corresponding arrows) indicate that 
left-right-left oscillation modes are persistent everywhere above the 
gi2 line, while in-out-in modes are persistent only above the higher 
^12 line. The spatially symmetric instability (^f 2 line) is the one rele- 
vant for experimental situations with symmetric traps. Squares mark 
values used in the dynamical simulations of Figs. [3] and IH (Table 1). 



modes describe in-phase motion of the two components and 
simply correspond to the excitation spectrum of a single- 
component BEC with interaction constant gu -\- gi2- Odd 
modes are more interesting — they describe out-of-phase mo- 
tion of two components and are therefore reflected in the spin 
dynamics. Additionally, due to the spatial inversion symme- 
try x —X, the solutions will also have well-defined spatial 
parity, and we can distinguish spatially symmetric and anti- 
symmetric modes. 

Typical eigenspectra are presented in Fig. [T] In the case 
of a tight trap A = 0.2, we notice two modes whose fre- 
quencies are nearly constant. These are even modes encod- 
ing single-component or in-phase physics. The lower one is 
the dipole (Kohn) mode with frequency equal to the trap fre- 
quency A. The second nearly-constant mode is the breath- 
ing mode, which for elongated traps takes value close to 
cj^ = 3A^. The breathing mode (oscillations of cloud size) 
is visible in the plots of Fig. [3l (Section [Tvb as a fast oscilla- 
tion of the total condensate widths. 

The two lowest-lying eigenmodes are odd modes encod- 
ing out-of-phase physics. For gi2 > 1, their frequencies are 
significantly below the breathing mode, and therefore lead to 
relatively slow oscillations in the spin density. This will also 
be visible in the real-time dynamics presented in Section HVl 
(first two columns of Figs. [3]and|4l). The forms of the corre- 
sponding eigenvectors are shown in the lower right of Fig.[T] 
The nature of the eigenvectors shows that the motion related 
to the lowest mode corresponds to the left-right oscillations of 
the two species, while the next odd mode corresponds to spa- 
tially symmetric spin motion. The frequencies of these two 
modes become imaginary at certain values of ^12, thus lead- 
ing to the onset of modulation instabilities. The antisymmetric 
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mode becomes unstable at smaller value of gi2 (gi2 ^1.6 for 
A = 0.2) in comparison to the symmetric mode (gl2 ^2.4 
for A = 0.2). In a spatially symmetric trap, there is no natural 
mechanism for exciting the spatially antisymmetric mode. On 
the other hand, any difference between gn and ^22 naturally 
excites the second (spatially symmetric) mode. Thus, the sec- 
ond mode, occurring at larger ^12, is the relevant instability 
for understanding the dynamics observed in experiments and 
explored numerically in Section Hvl 

We find similar excitation spectra for trap strengths A span- 
ning several orders of magnitude. The spatially antisymmetric 
mode becomes unstable before the spatially symmetric mode, 
and both instabilities get closer to 1 as the trap gets shallower. 
For example, for A = 10~^ (also shown in Fig.[T]) the lowest 
instability sets in for ^ 1-02, while the next one appears 

9 12 ^ ^'Ob. The distinction between two instabilities be- 
comes ever smaller as we go toward a uniform system A ^ 0, 
where the phase- separation condition Eq. ([T]) becomes exact. 
Nevertheless, even for shallow traps, the issue is not purely 
academic as the precision in experimental measurement and 
control of scattering lengths continues to improve 10, H^. 

In Fig. [2] (main panel), the results of the stability analyses 
are combined to present a dynamical "phase diagram". The 
two lines show the two instabilities (^^2 9 12) ^ function 
of trap strength A. For very shallow traps, the two transition 
lines merge as gl2 ^ 9i2 ^ 1- The lower transition line (^^2) 
was previously introduced in Ref . |[T4l] . However, for a trap 
and initial state with left-right spatial symmetry, this is not the 
relevant dynamical transition line, because the first even mode 
only becomes unstable at some higher gi2 value, given by the 
9i2 line. 

In the next Section, we will see that spin pattern dynamics 
is indeed only generated when the inter-component repulsion 
gi2 exceeds the second instability line (gi2 > ^12), and that 
crossing the first instability (^^2 < ^12 < ^12) enough 
for pattern formation in a spatially symmetric trap. 

IV. DYNAMICAL FEATURES ACROSS THE PARAMETER 
SPACE 

In this Section we present and analyze the dynamics 
obtained from direct numerical simulation of the Gross- 
Pitaevskii equations (0 and ©, after the system is initially 
prepared in the ground state of the situation gn = ^22 = 
gi2 = 1. The subsequent dynamics is performed with gn = 
1, 922 = 1-01, and several different values of gi2 for each trap 
strength A. 

It is difficult to show the full richness of pattern dynam- 
ics through plots of a few quantities. We choose to show the 
dynamics through two types of plots (Figs. [3] and |4]). Fig. [3] 
shows the time dependence of the root mean square widths of 
the two components 

/oo 
x^\i;i,2{x,t)\''dx, (6) 
-00 

while Fig.|4]shows density plots of the density difference (spin 
density), \ipi{x^t)\'^ — \ilj2{xjt)\'^ . In both figures, each row 
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TABLE L Parameters from the first five columns are used for the 
plots in Figs.[3]and Figs.lH The instability values gi2 and gl2 (intro- 
duced in Figs. [T]and[2l and discussed in Section|lIll) are also given for 
each trap strength. 



corresponds to a different trap strength (A), and we approach 
the shallow trap (Thomas-Fermi) limit going from top to bot- 
tom. 

For each A the four values of gi2 from Table [Bare used for 
Figs. [3] and m We have chosen gi2 values such that the first 
panel in each row is in the parameter region where there are 
no instabilities, the second one is in the region where the only 
instability is the antisymmetric one, and the third on each row 
is at gi2 values just above the second, relevant, instability. The 
fourth panel on each row is at higher gi2 values. The choice 
of gi2 values with respect to instability lines is clear in the 
tighter traps of the top three rows, as also shown by squares in 
Fig. [21 For shallow traps (lower rows), the instability lines are 
too close together and too close to gi2 = 1, so making such 
choices is not meaningful. In the following, as we compare 
features of the different columns, we implicitly exclude the 
lowest row (smallest A). This is also indicated by the fact that 
the schematic instability lines in Figs.[3]and[4lare not extended 
to the lowest row. 

Broadly speaking, we note that there is only regular 
(collective-mode) dynamics in the second-column figures 
(^12 < ^12 < ^12) ^^^^ though an instability is present. There 
is generally a sharp difference between the second and third 
figure in each row, indicating that the second instability {gl2) 
is the relevant one. The fourth panel on each row is at higher 
gi2 values, showing more rich dynamics. 

In Fig. [31 we show time-dependence of the individual 
widths {wi, W2) and also of the total root mean square width, 
w{t) = ^ i^i^iit) + W2{t))/2. Consistent with our observa- 
tion that spatially symmetric modes (and not the antisymmet- 
ric ones) are naturally excited in the current setup, the dy- 
namics shows signatures of the two most prominent spatially 
symmetric modes noted in Fig. [T] The breathing mode is the 
easiest to notice and most ubiquitous — it shows up in almost 
every parameter choice as oscillations in the total density (in- 
phase in the two components), with a typical period given by 
27r/>/3A ^ 3. 63/ A. This follows from the frequency of this 
mode being almost constant near ^/3X. 

We also see out-of-phase motion of the two components, 
associated with the lower spatially- symmetric mode in Fig. 
[TJ which has odd species parity. In the first two columns of 
Fig. [3 corresponding to smaller values of ^12 such that this 
mode has small real frequencies, this is excited as a regular 
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FIG. 3. (Color online.) Time evolution of root-mean- square widths after 7r/2 pulse (interaction quench). First component width wi{t) 
is shown as blue dashed line, second component width W2{t) is shown as red solid line (gray solid without color), the total width w(t) — 
^/ (w\(t) + i(;|(t))/2 is the black solid line intermediate between the other two. From top to bottom: tight to shallow traps. For each trap 
strength, four values of gi2 (indicated near top of each panel) from Table IJ are used. The two lines separating first and second column (red 
dashed) and second and third column (black solid) indicate the 'positions' of instability lines, from Figure|2] While the first two columns look 
qualitatively the same and show regular oscillatory dynamics, in the third column we observe aperiodic motion of stronger amplitude that we 
relate to the onset of spin pattern dynamics. The spin dynamics is even more pronounced in the fourth column. 



'spin' mode. For example, at A = 10~^ and gi2 = 1.04, 
we observe an out-of-phase oscillation with the period of ap- 
proximately ~ 30, much slower than the breathing mode. In 
addition, the oscillation period of the out-of-phase motion is 
slower in the second than in the first column of each row, cor- 
responding to the decreasing frequency of the mode, as seen 
in stability analysis (Fig.[T]). Once gi2 becomes large enough 
that the instability threshold for this mode is crossed, the os- 
cillation amplitudes increase sharply and the width dynam- 
ics becomes strongly aperiodic and irregular (third column of 
Fig. O. This signifies the onset of pattern dynamics, as op- 
posed to the excitation of a regular collective mode around a 
stable state. Irregularity of the width dynamics at stronger gi2 
is even more apparent in the fourth column of Fig. [3] 

It is noteworthy that the spatially antisymmetric modes play 
no role and do not show up in these dynamical simulations. 
We see no signature of the Kohn mode. Nor do we see any 
sharp change associated with the instability of the antisym- 
metric mode, i.e., there is no sharp difference between the first 
two columns of Fig. [S] 

In Fig. [4] we show the dynamics of the "spin density" 



— |?/^2(^,^)p. The case of very shallow traps (last 
row), resembles the data in Refs. ifTvl fTsll. As in Fig.[3l the 
first two columns show regular oscillations, corresponding 
to collective modes without instability. A sharp change oc- 
curs, not across the first instability line (between 1st and 2nd 
column), but instead across the second instability line (2nd 
and 3rd columns), especially for tighter traps (top three rows) 
where comparison with instability lines is meaningful. The 
sharp change can be noted through the color scales, which is 
dramatically different between second and third columns in 
the upper rows. 



V. LENGTH SCALES OF PATTERNS 

In homogeneous stability analysis, the length scale of pat- 
terns is inferred from the wavevector (momentum) at which an 
instability first occurs. Since we perform our stability analysis 
specifically for trapped systems, we do not have a momentum 
quantum number. Nevertheless, the eigenvectors of the un- 
stable modes contain information about the form of patterns 
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FIG. 4. (Color online.) Spin dynamics subsequent to the 7r/2 protocol, shown via the density difference A (|^i(x, t)|^ — |'02(x, t)|^). 
Traps and gi2 values are the same as in Fig. [3] and Table HI A decreases from 10~^ to 10~^ from top to bottom and gi2 values are indicated 
near top of each panel. As in Fig. [3] the black solid line and the red dashed line indicate the instability lines from the "phase diagram" of 
Fig. 12] Note the sharp change of color-scale ranges between second and third columns in the upper rows, indicating that the dynamics changes 
dramatically only across the second instability line. 



generated in the dynamics of the trapped system. This is il- 
lustrated in Fig. [51 where the eigenvectors of the lowest unsta- 
ble even-parity modes are shown for several values of gi2, to- 
gether with the spin patterns generated in the non-equilibrium 
dynamics. There is a close match between the distance be- 
tween nodes of the eigenvectors (rough analog of 'wavevec- 
tor') and the length scales involved in the patterns. 



In Fig.|4]we see that the patterns contain more spatial struc- 
ture in shallow traps. The top two rows (tight traps) only 
show in-out-in type of patterns. This can be understood from 
the idea that the interactions induce length scales ('healing 
lengths') in the problem, which are smaller for larger interac- 
tions, and which set the length scale of spatial structures. For 
tight traps, the healing length set by the interactions is large 
or comparable to the cloud size, so that only global dynamical 
patterns are generated. In such traps, generation of complex 
patterns with many spatial oscillations would require much 
higher values of (^12/^11 — 1). For shallow traps, the healing 
length becomes much smaller than the cloud size; as a result 
one can have a multitude of dynamical spin structures in the 
system, of the type seen in experiments and prior simulations 
iSSllSl]- This heuristic explanation can be made quantitative 
by counting the number of nodes appearing in the eigenmodes 
(as in Fig.O. 



VI. CONCLUSIONS AND OPEN PROBLEMS 

In this article we have analyzed a widely used dynamical 
protocol for two-component BECs, which involves starting 
from the ground state of one component and switching half 
the atoms to a different component through a 7r/2 pulse. We 
have presented a stability analysis suitable to the trapped sit- 
uation, and also presented results from extensive dynamical 
simulations. Through an analysis of unstable modes, we have 
presented a classification of the parameter space into a number 
of dynamically distinct regions, in relation to the prototypical 
initial state. This may be regarded as a dynamical "phase dia- 
gram". 

In the 'stable' regime of parameter space (no modulation 
instabilities), our stability analysis explains the observed slow 
spin oscillations compared to the fast breathing mode oscilla- 
tions of the total density. We demonstrate that the important 
"phase transition" line for spatially symmetric situations rele- 
vant to most experiments is not the first instability (studied in 
Ref. 1IT4I1 ). but a second transition line. The first instability is 
antisymmetric in space, and as a result is not naturally excited 
in a symmetric trap. 

Our stability analysis is performed relative to a stationary 
state of the situation gu = ^22- The 7r/2 pulse of the ex- 
periments (in the cases where ^117^^22) can be considered as 
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FIG. 5. (Color online.) Top: Eigenvectors of the most unstable 
even eigenmodes, from the stability analysis of Section [nil for A = 
10"^, and gi2 = 1.02, 1.06, and 1.12 from left to right. Below each 
eigenvector, the corresponding spin dynamics after the 7r/2 protocol 
(parameters of Section HVl) is shown through the time evolution of 

|V^l(x,t)|2-|V^2(x,t)|^ 



turning on a nonzero (^n — ^22), i.e., turning on 'buoyancy' 
such that one component gains more energy by being in the 
interior of the trap compared to the other. This helps to select 
instability modes which are symmetric in space. 

Since we have used a stability analysis with gn = ^22 to 
analyze dynamics with ^117^^22, an obvious question is how 
the ratio ^22/^11 affects the regime of applicability of this 
scheme. We expect that features of this (gu = ^22) stabil- 
ity analysis are useful for dynamical predictions as long as 
^12/^11 — 1 is roughly more than ^22/^11 — !. For example, for 
shallow traps (small A), the instabilities occur at ^12/^11 — 1 
values comparable to 0.01, which is why the placement of pa- 
rameters in the three dynamical regions of the 'phase diagram' 
(Fig. O is not meaningful for the smallest A values (lowest 
rows ofFigs.[3]and|4]). 

For the stability analysis we used a reference stationary 
state which is of course not the initial state: the initial state is 
the ground state for gu = ^22 = 912, while the reference state 
is the lowest-energy spatially symmetric stationary state cor- 
responding to the final value of gi2. The instability lines found 
in this stability analysis would describe even better a situation 
where the dynamics is triggered by a small quench of ^12, as 
opposed to the changes of gi2 that we consider here, which 
can be relatively large. We have looked at some examples of 
this type of dynamics and indeed find instabilities matching 
the stability analysis extremely well. However, although the 



initial state in the 7r/2 dynamics is somewhat different from 
the reference state of our stability analysis, our results show 
that this stability analysis does provide an excellent overall 
picture of the dynamics generated by the 7r/2 protocol. 

Our work opens up a number of questions deserving further 
study. First of all, we have thoroughly explored the A-^12 pa- 
rameter space, while assuming that the intra-component inter- 
actions gu and ^22 are unequal but close in value. The regime 
of large difference (^n — ^22) clearly might have other inter- 
esting dynamical features which are yet to be explored. 

Second, in this work we have restricted ourselves to the 
mean field regime. While the mean field description cap- 
tures well the richness of pattern formation phenomena (c.f. 
Refs. HHIitIIiIIIIS] in addition to present work), it may 
be worth asking whether quantum effects beyond mean field 
might have interesting consequences for the patter dynam- 
ics generated by a 7r/2 pulse. For bosons in elongated traps, 
regimes other than mean field (such as Lieb-Liniger or Tonks 
regimes) may occur naturally in experiments [37-41]. Dy- 
namics subsequent to a 7r/2 pulse in strongly interacting ID 
gases outside the mean field regime is an open area of investi- 
gation. 

Third, we have assumed a spatially symmetric trap and an 
initial condition with spatial symmetry, and this plays a cru- 
cial role in the selection of instability channels. In a real- 
life experiment, the trap will have some left-right asymmetry. 
Also, thermal and quantum fluctuations can initiate spatially 
antisymmetric excitations. The extent to which a small spatial 
asymmetry affects spin dynamics remains unexplored; in such 
a case we would have som e typ e of competition between two 
types of instabilities. Ref. (]M has studied dynamical effects 
of fluctuations (noise), but the effects of thermal and quan- 
tum fluctuations is yet to be studied in the context of a 7r/2 
protocol. 

Finally, one could consider time evolution and spatiotem- 
poral patterns generated by a 7r/2 pulse in the presence of an 
optical lattice, described by the dynamics of a two-component 
Bose-Hubbard model. This is a situation easy to imagine real- 
izing experimentally. One could speculate complex interplay 
between spin dynamics and the spatial arrangement of Mott 
and superfluid regions. 
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